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Abstract 

The Wightman function and the vacuum expectation values of the held squared and 
of the energy-momentum tensor are obtained, for a massive scalar held with an arbitrary 
curvature coupling parameter, in the region between two inhnite parallel plates, on the 
background of de Sitter spacetime. The held is prepared in the Bunch-Davies vacuum state 
and is constrained to satisfy Robin boundary conditions on the plates. For the calculation, a 
mode-summation method is used, supplemented with a variant of the generalized Abel-Plana 
formula. This allows to explicitly extract the contributions to the expectation values which 
come from each single boundary, and to expand the second-plate-induced part in terms 
of exponentially convergent integrals. Several limiting cases of interest are then studied. 
Moreover, the Casimir forces acting on the plates are evaluated, and it is shown that the 
curvature of the background spacetime decisively influences the behavior of these forces at 
separations larger than the curvature scale of de Sitter spacetime. In terms of the curvature 
coupling parameter and the mass of the held, two very different regimes are realized, which 
exhibit monotonic and oscillatory behavior of the vacuum expectation values, respectively. 
The decay of the Casimir force at large plate separation is shown to be power-law (monotonic 
or oscillating), with independence of the value of the held mass. 



PACS numbers: 04.62. +v, 04.20. Gz, 04.50.-h, ll.10.Kk 



1 Introduction 

The Casimir effect [1] is now known to be common to systems of very different kind, involving 
fluctuating quantities on which external boundary conditions are imposed. It can have important 
implications on all scales, from subnuclear to cosmological. Imposing boundary conditions on a 
quantum field leads to a modification of the spectrum of zero-point fluctuations and results in 
the shifting in the vacuum expectation values for physical quantities, such as the energy density 
and stresses. In particular, the confinement of quantum fluctuations induces forces that act on 
the constraining boundaries. The particular features of the resulting vacuum forces depend on 
the nature of the quantum field, on the type of the spacetime manifold, the geometry of the 
boundaries, and on the specific boundary conditions imposed on the field. 
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An interesting topic in the investigation of the Casimir effect is its explicit dependence on 
the geometry of the background spacetime. As usual, the relevant information is encoded in 
the vacuum fluctuations spectrum and, not surprisingly, analytic solutions can be found for 
highly symmetric geometries only. In special, motivated by Randall-Sundrum type braneworld 
scenarios, investigations of the Casimir effect in anti-de Sitter (AdS) spacetime have attracted 
a great deal of attention. The braneworld corresponds to a manifold with boundaries and all 
fields which propagate in the bulk will give Casimir-type contributions to the vacuum energy 
and, as a result, to the vacuum forces acting on the branes. The Casimir effect provides in this 
context a natural mechanism for stabilizing the radion field in the Randall-Sundrum model, as 
required for a complete solution of the hierarchy problem. In addition, the Casimir energy gives 
a contribution to both the brane and bulk cosmological constants and, hence, it has to be taken 
into account in any self-consistent formulation of the braneworld dynamics. The Casimir energy 
and corresponding Casimir forces for two parallel branes in AdS spacetime have been evaluated 
in Refs. [2] by using either dimensional or zeta function regularization methods. Local Casimir 
densities were considered in Refs. pUS]. The Casimir effect in higher-dimensional generalizations 
of the AdS spacetime with compact internal spaces has been investigated in [5]. 

Another popular background in gravitational physics is de Sitter (dS) spacetime. Quantum 
field theory in this background has been extensively studied during the past two decades. Much 
of the early interest was motivated by questions related with the quantization of fields propa- 
gating on curved backgrounds. The dS spacetime has a high degree of symmetry and numerous 
physical problems are exactly solvable on this background. Importance of such theoretical work 
was increased with the appearance of the inflationary cosmology scenario [6]. In most infla- 
tionary models, an approximate dS spacetime is employed with the aim to solve a number of 
problems in standard cosmology. During the inflationary epoch, quantum fluctuations in the 
inflaton field introduce inhomogeneities which play a central role in the generation of cosmic 
structures from inflation. More recently, astronomical observations of standard-candle super- 
novae, galaxy clusters and the cosmic microwave background [7] have clearly indicated that 
at present our (local) universe is accelerating and can be well approximated by ACDM, FRW 
cosmology with a positive cosmological constant A. Now, if the universe, as it seems, is going 
to accelerate for ever, this cosmology will lead asymptotically to a dS universe. Another moti- 
vation for the investigation of dS-based quantum theories is related to the holographic duality 
known to hold between quantum gravity on dS spacetime and a quantum field theory living on 
its boundary, identified with the timelike infinity surface of the dS spacetime. 

Motivated by the above considerations — and since they are ingredients that more full-fledged 
models will necessarily have to incorporate — we will here calculate the Casimir densities and 
forces arising for the geometry of two parallel plates on the background of (D + l)-dimensional 
dS spacetime. Previously, the Casimir effect on the background of dS spacetime described in 
planar coordinates was investigated in Refs. [8] for a conformally coupled massless scalar field. 
In this last case the problem is conformally related to the corresponding problem in Minkowski 
spacetime and the vacuum characteristics are generated from those for the Minkowski counter- 
part, just by multiplying with the conformal factor. In particular, for the geometry of a single 
plate, the vacuum expectation value of the energy-momentum tensor vanishes. The Casimir 
density induced by a single plate for a massive scalar field with an arbitrary curvature coupling 
parameter has been considered in [9]. In [10] the vacuum expectation value of the energy- 
momentum tensor for a conformally coupled scalar field was investigated in dS spacetime with 
static coordinates in presence of curved branes, on which the field obeys the Robin boundary 
conditions with coordinate dependent coefficients. In those papers the conformal relation be- 
tween dS and Rindler spacetimes and the results for the Rindler counterpart were used. More 
recently, the Casimir density in a dS spacetime with toroidally compactified spatial dimensions 
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has been investigated in 

The outline of the paper is as follows. In the next section the positive frequency Wightman 
function will be evaluated for a scalar field with general curvature coupling parameter with 
Robin boundary conditions on two parallel plane boundaries in the background of dS spacetime. 
Among the most important quantities describing the local properties of a quantum field and 
the corresponding quantum back-reaction effects are the expectation values of the field squared 
and of the energy-momentum tensor. These quantities, in the region between the plates, will be 
investigated in Sects. [3] and HI Simple asymptotic formulae are obtained both for small and for 
large plate separations. The Casimir forces acting on the plates are obtained in Sect.EJ Finally, 
Sect. [6] contains a summary of the work done together with an outlook. 

2 The Wightman function in de Sitter spacetime with two par- 
allel plates 

Consider a quantum scalar field (f{x) on a (D + l)-dimensional dS spacetime background, as 
coming from a cosmological theory, as described in Sect. 1, with positive cosmological constant, 
A. The corresponding line element can be written in planar (inflationary) coordinates, which 
are most appropriate for cosmological applications: 

ds 2 = dt 2 -e 2t / a J2( dzi ) 2 - (!) 
i=i 

Here, the parameter a is related to the cosmological constant through the expression a = 
D(D — 1)/(2A). For the discussion to follow, in addition to the synchronous time coordinate, 
t, it is very convenient to introduce the conformal time, r, defined as r = -ae~^ Q , — oo < 
t < 0. In terms of this coordinate the metric tensor takes the conformally flat form: = 
(a/r) 2 diag(l,-l,...,-l). 

The dynamics of a massive scalar field with curvature coupling parameter are governed by 
the Klein-Gordon equation 

(v t V l + m 2 + (R^ip = 0, (2) 

where V; is the covariant derivative operator and R = D(D + l)/a 2 is the Ricci scalar for dS 
spacetime. The special cases £ = and £ = £d = (D — 1)/4D correspond to minimally and to 
conformally coupled fields, respectively. The importance of these two special cases comes from 
the fact that, in the massless limit, the corresponding fields mimic the behavior of gravitons 
and photons, respectively. Note that non-minimal coupling is required by the renormalizability 
condition for interacting theories in curved spacetime |12| . 

In this paper we will be interested in the study of the Casimir densities and of the mutual 
forces occurring for the geometry of two infinite, parallel plates in dS spacetime. The plates are 
located at z D = a\ and z D = a2, a\ < a^. As the most general set up, we assume that on these 
boundaries the scalar field obeys Robin boundary conditions (BCs) 

(l + P j n l ViMx) = [l + p j (-iy- 1 d D } V (x)=0, z D = aj , j = 1,2, (3) 

with constant coefficients f3j. For /3j = these BCs reduce to Dirichlet ones, and for /3j = oo 
to Neumann BCs. The choice of different boundary conditions on the plates corresponds, in 
physical terms, to using different materials for the same. The imposition of BCs leads to a 
modification of the vacuum expectation values (VEVs) for physical quantities, as compared 
with those in the situation without boundaries. In the discussion below we will assume that the 
quantum scalar field is prepared in a dS invariant Bunch-Davies vacuum state |13j . 
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Among the most important characteristics of the vacuum state are the VEVs of the field 
squared and of the energy-momentum tensor. These VEVs are obtained from the corresponding 
positive frequency Wightman function W(x,x') in the coincidence limit of the arguments. The 
Wightman function is also of the essence for the consideration of the response of particle detectors 
at a given state of motion (see, for instance, [H]). Expanding the field operator over a complete 
set {(p a (x), (p^(x)} of solutions to the classical field equation, satisfying the boundary conditions, 
the positive frequency Wightman function is best expressed as the mode-sum 

W(x,x') = {0\<p(xMx')\0) = ^tp.Wtitf), (4) 

17 

where the collective index a labels the solutions. 

In the region between the plates, a\ < z D < a 2 , the eigenfunctions realizing the Bunch- 
Davies vacuum state and satisfying the BC on the plate at z D = a\, have the form 



a (x) = C aV D / 2 Hi 1 \ V K) C os[k D (z D -a^ + a^ko)]^, (5) 



with the notations 77 = |r|, K = \/k 2 + k^, and 

''/3ix - 1 , . 

InEq. ©, z= (z 1 , . . . , z - 1 ) is the position vector along the dimensions parallel to the plates, 
k = . . . , kz>-i), and the order of the Hankel function H„ (x) is given by 

v = [D 2 /4 - D(D + 1)^ - m 2 a 2 } 1/2 . (7) 

Note that v is either real and nonnegative or purely imaginary. For a conformally coupled 
massless field v = 1/2 and the function H„\x) in ([5|) is expressed in terms of elementary 
functions. From the boundary condition on the second plate z D = a? we find that the eigenvalues 
for kr> are solutions of the transcendental equation 

(1 - hb 2 y 2 ) sin y- (pi + b 2 )y cos y = 0, (8) 

where y = k^a and bj = fij/a, with a = a 2 — a\ being the separation between the plates. In 
the discussion below we will assume that all zeros are real. In particular, this is the case for 
the conditions bj ^ (see [E]). The positive solutions of Eq. (Jl]) will be denoted by y = A n , 
n = 1, 2, . . ., and for the eigenvalues of kn one has ko = X n /a. The eigenfunctions are specified 
by a set of eigenfunctions a = (k, n). 

The coefficient C„ in ([5]) is determined from the orthonormalization condition 

/ 1-0,2 
dz / dz D ^/\g~\9™ [<P*{x)d T y*Ax) - yl,{x)d T y a {x)] = i5 nn ,5(k - k'). (9) 
J ai 

By using the Wronskian relation for the Hankel functions, we find 

C 2 = (2n) 2 ~°gi- D e^-^/ 2 
u 4a{l + cos[A n + 2ai(A n /a)]sin(A n )/A n }' 1 ' 

the star meaning complex conjugate. 
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Combining Eqs. (HD, ©, (fTUj) . for the Wightman function in the region between the plates, 
we find 



W(x,x') 



4(r]r]') D / 2 
{2ir) D aa D - 

cos[X n (z D 



/oo 
dkel k.(z-z') J2 K„(T]k n e~ 7ri/2 )K„(r]' k n e ni/2 ) 
n=l 



ai)/a + ai(A n /a)] cos[\ n {z Dl - ai)/a + ai(\ n /a)] 
1 + cos[A n + 2«i(A n /a)] sin(A n )/A n 



(11) 



with k n = \Jk 2 + A 2 /a 2 and where we have written the Hankel functions in terms of the modified 
Bessel function K v {x). It is well known that in dS spacetime without boundaries the Bunch- 
Davies vacuum state is not a physically realizable state for Rev ^ D/2. The corresponding 
Wightman function contains infrared divergences arising from long- wavelength modes. In the 
presence of boundaries, the BCs on the quantized field may exclude these modes and the Bunch- 
Davies vacuum becomes a realizable state. An example of this type of situation is provided by 
the geometry of two parallel plates described above. In the region between the plates and for 
BCs with f3j ^ 0, j3j ^ oo, there is a maximum wavelength, 27ro/Ai, and the two-point function 
([lip contains no infrared divergences. Mathematically, this situation corresponds to the one 
where in the arguments of the modified Bessel functions we have k n ^ \\/a. 

As we do not know the explicit expression for A n as a function of n, and being the summand 
in (jlip a strongly oscillating function for large values of n, this formula is not convenient for the 
evaluation of the VEVs of the field squared and of the energy-momentum tensor. In order to 
obtain a useful alternative representation, we apply to the series on n the summation formula 

dais] 



7rA n /(A r 



7T 



/(0) 



^ A n + sin(A n ) cos[A n + 2a 1 (X n /a)] 



2 1 



+ 



+i 



dz 



h 

f(iz 



dzf(z) 



(h*-i)(fa*-i) - 2g 



1 



(12) 



being 



K v (j}e-™l 2 ^k 2 + z 2 /a 2 )K u (r/e™ /2 ^k 2 + z 2 /a 2 ) 

x cos[z(z D — ai)/a + a\(z/a)] cos[z(z Dl — a{)/a + a\{z/a)\. 



(13) 



In the case bj > this function has poles on the imaginary axis and the corresponding residue 
terms should be added to the right hand side of the summation formula. In order to easy 
the presentation, in the discussion below we will just consider the case bj ^ 0, but a similar 
procedure is valid in the general case. 

Use of the summation formula (|12p with (|13p allows us to write the Wightman function in 
the decomposed form 

W(x,x') 



W l (x,x') + 



2a 



l-D 



du 



(2tt) d 
cosh[u(z D 



dkei k-(z-z') 

di) + <5i(it)] cosh[u(z D ' — ai) + ai(w)] 



xy 



-D 



ci(u)C2 (u)e 2au 
I v {r\ y)K v {r]y) + l v {j]y)k v {r\ y) 



1 



(14) 



where the function ol\{u) is defined by the relation e 2ai ^ 
notations 

KM = y D ' 2 K v (y), Uv) = V° 12 Mv) + I-v{v)\ , 



y=Vu 2 ~k 2 

ci(u) and we have introduced the 



(15) 
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and 



Note that one has Cj(u) 
Eq. ([IT 



Cj(u) 



/3 jU - 1 
JIu + 1 ' 



(16) 



-1 in the case of Dirichlet BCs and Cj(u) = 1 for Neumann BCs. In 



Wi(x,a/) 



8( W ') D / 2 
{2ix) D + 1 a D - 1 

x cos[u{z 



dke ik-(z-z') 



cZu cos[u(z Z) — a\) + ai(u)] 



oi) + 



ai(u)]^(? ? v / A; 2 + !i 2 e-" /2 )^(r,Vit 2 + ^e™/ 2 ) (17) 



is the Wightman function corresponding to a single plate at z = a\. 

The VEVs for the geometry of a single plate were investigated in [9j . Denoting the Wightman 
function for the dS spacetime without boundaries by Wds(x,x'), the corresponding Wightman 
function for a plate located at z D = aj can be written in the form 

W d (x,x') = W dS (x,x') + W<. l) (x,x'), 
where the part induced by the plate is given by the expression 



(18) 



wP(x,x') 



a 



l-D 



2(2tt) d 

-D 



dke4 k-(z-z') 



du • 



-u\z D +z D '-2a :j \ 
Cj(u) 



xy 



I v (rj 'y)k y {r\y) + I v {ryy)k v {r] y) 



(19) 



y=V u' 2 — k 2 

The two-point function in the dS spacetime without boundaries was investigated in [HI [18] (see 
). It is given by the formula 

» l - D T(D/2 + v)T{D/2-i 



also 



Wds(x,x') 



2 (0+3)/2 7r (D+l)/2 ( n 2 _ I 



,(-D-l)/4 J iz-1/2 



D (l-£>)/2/ \ 



where Pv{x) is the associated Legendre function of the first kind and 



u 



-1 + 



-*") 2 - to -*/) 

2r/r/' 



/\2 



(20) 



(21) 



By using the definition of cti(u), the Wightman function in the region between the plates 
can be written in the more symmetric form: 

..l-D 



W(x,x') 



W dS (x,x')+ Wf{x,x') + 



V 

i=i,2 



o 



00 2cosh[u(z D 



2(2vr) D 

+ Ej = l,2 



dke ik-(z-z') 



-u\z L 



-2a j | 



/ c i(«) 



xy 



-D 



ci(u)c2(u)e 2au 
I v {r\ y)k v {r\y) + I u (r)y)K u (r)'y) 



1 



(22) 



i/=\/ u 2 — k 2 

where the last term on the rhs can be referred to as the interference part. This integral repre- 
sentation of the Wightman function is valid for Rei^ < 1. As it has been shown in Ref. |17j . 
the quantized graviton field in D = 3 dS spacetime is equivalent to a pair of minimally coupled 
massless scalar fields. Thus the representation (|22p does not apply to the graviton. In the region 
z D < a\ (z D > 02) the Wightman function coincides with the corresponding function for a single 



plate located at z 



D 



ai [z 



D 



0,2) and is given by the expression ^8[) . with j = 1 (j = 2). 



The results obtained in the present paper can be applied to a more general problem where the 
cosmological constant is different in separate regions z D < ai, a\ < z D < 0,2, and z D > 0,2- In 
this case the plate can be considered as a simple model of a thin domain wall separating the 
regions with different dS vacua. 
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3 Vacuum expectation value of the field squared 



Once we have the Wightman function, we can proceed to evaluate the VEV of the field squared, 
by taking the coincidence limit of the arguments. In this limit the Wightman function is di- 
vergent and some renormalization procedure is needed. The important point here is that for 
points far away from the boundaries the divergences are the same as for the dS spacetime with- 
out boundaries. As in our previous procedure we have already extracted, from the Wightman 
function, the part Wds{x,x'), the renormalization of the VEVs is just reduced to the renormal- 
ization of the part corresponding to the geometry without boundaries, which is already done in 
literature. For the further discussion of the VEVs in the coincidence limit it is convenient to 
use the notations 

( x 2 _ y 2)(D-3)/2 

Ci{x/r))c2{x/ri)e' 2ax l r i — l' 
g{P jU ,yu) = c J (u)e 2u y + e- 2u y/c J (u) + 2. (24) 

Introducing also the function 

F v {y) = y D Mv) + i-u{y)\ K u (y), (25) 

The VEV of the field squared can be expressed in the decomposed form 

/ 2\ / A , / 2\ A D f°° F u (y) 
{<p) = (<P >ds + (<p )j + -7531 / dy-o-T 

a Jo y 

x dx H(x,y)g(/3jx/rj,\z D - aj\x/r]), (26) 
Jy 

where 



4(4vr)-( D+1 )/ 2 

Ad = r((D-i)/2)' (27) 

In Eq. (]26p . (y 2 )dS is the renormalized VEV in dS spacetime without boundaries, and the part 

2 V A D r^Mv) ' X 



y 

D 



x ( x 2 _ y^-S)/2 e -2 X \z--a j \/ V/cj{x/v)f (28) 



ratio \z D -a j\/rj in the formulae for the VEV is the proper distance of the observation point 
from the plate at z D = aj measured in units of the curvature scale, a (note that a is the dS 



is induced by a single boundary located at Note that, due to the dS invariance of 

the Bunch-Davies vacuum state, the VEV (y? 2 )ds does not depend on the spacetime point. The 

horizon size). The boundary induced VEV is a function of the combinations \z D — a,j\/r] and 
fij/t] only. This property follows from the maximal symmetry of the dS spacetime. 

We can write the expression for the VEV of the field squared in a more symmetric form by 
using (122]) 

<^> = <Ais+ ^(A' + A(^ 2 >, (29) 
where the interference term is given by the expression 

A POO f'OO 

A ^ 2 ) = ^D R iJ Q dyy 1 - D F u (y)j dx H(x,y)h(x/ V , z D ), (30) 
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with the notation 

h(x, z D ) = 2+Y, e~ 2x \ zD - a ^ / Cj (x). (31) 
i=i,2 

Note that the interference part is finite everywhere, including the points on the plates. The 
surface divergences in the VEV are contained in the single plate parts only. For points near 
the plates the total VEV is dominated by these contributions. In particular, near the plate at 
z D = a,j one has ((p 2 ) ~ (\z D — aj\/rf) l ~ D . The corresponding result for two parallel plates in 
Minkowski spacetime is obtained from (|30p in the limit a — > oo. In this limit one has v ~ ima 
and the modulus of the order of the modified Bessel function is large. In addition, we have 
77 — > a. Details of the corresponding limiting transition are given in Appendix, for the explicit 
case of the vacuum forces acting on the plates. 

For a conformally coupled massless scalar field (£ = £d, m = 0) one has v = 1/2 and 
F v {y) = y D ~ l . In this case, for the interference term we find 

{iP ' ~ (4vr)^/2r(D/2) J Cl (x)c 2 (x)e^ - 1 ' ^ > 

This result (1321) could also have been directly obtained by using the fact that, in the special case 
under consideration, the problem is conformally related to the corresponding one for Robin plates 
in Minkowski spacetime [151 H6j. From this relation, it follows that A(</? 2 ) = a l ~ D (r])A{(p 2 )( M \ 
with scale factor a{rj) = a/rj, which leads to the result (|32p . 

Formula (130|) is further simplified for Dirichlet and Neumann BCs. Using the expansion 
(e 2 — 1) _1 = Yln=i e ~ nz an d explicitly integrating over x, we find 



n = 

2/c/ a -i(2noy) - 5^ ^ f D /z-i(?V (™ + V° ~ %l)) 
3=1,2 



(33) 



where fu(x) = K^(x)/x^, J=D,N for Dirichlet and Neumann boundary conditions, 5^ = 1, 
6 {N) = _ L Note that 

as in the case of the Minkowski bulk, the coordinate dependent parts 
in the VEV have opposite signs for Dirichlet and Neumann scalars. By taking into account 
that for real values ^ v < 1 the function F u (rjy) is non-negative and the function fu{%) is 
monotonically decreasing, from (|33p we conclude that the interference term for Dirichlet and 
Neumann scalars is always positive. 

Now we turn to the investigation of the interference part in the VEV of the field squared 
in the asymptotic regions of the ratio a/77. This ratio is the proper distance between the plates 
in units of the dS curvature scale a. For small values of a/77 the main contribution to (f30j) 
comes from large values of y, for which we have F u (y) ~ y D ~ l ■ As a result, at leading order the 
interference part in the VEV of the field squared coincides with the corresponding quantity for 
a conformally coupled massless field, and is given by Eq. (f32j) . In particular, this part is positive 
for Dirichlet or Neumann BCs on both plates, and is negative for Dirichlet BC in one plate and 
Neumann on the other. 

For large proper distances between the plates one has a/r] ^> 1. In order to find the leading 
terms in the corresponding asymptotic expansion we introduce in (I30p new integration variables 
u = ax/r/ and v = ay/rj. In these variables the argument of the function (|25p becomes vq/a and 
in the limit under consideration it is small. By using the asymptotic formulae for the modified 
Bessel functions to the leading order, we have 



rn 



y « 1, (34) 
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where a v = 1 and a v = 2 for real and imaginary u, respectively. By taking into account this 
relation, as the next step we apply the integration formula 



00 dv yi ~ 2u r du {u2 - v2){D ~ 3)/2f{u) = r( 2 i rag ) +ij/2~- i i / ) 2) r dr rD ^ if{r) > (35) 



for a given function f(u). Assuming also that \/3j\/a <C 1, for the remaining integral we can use 
the expression 

00 x D-2u-l e -x\z D - aj \/a ^ » {5 1 5 2 ) n T{D -2v) 

«5i«5 2 e- - 1 " ^ ( n + |^_ a .|/ a )D-2,' ( 36 ) 

where <5j = c, (0). Note that 5j = — 1 for non-Neumann BC (|/3j| < 00) on the plate at ctj, while 
5j = 1 if the BC is Neumann. Finally, by using the duplication formula for the gamma function 
in ()36p . we find the following asymptotic behavior 

..1--D 



(4vr) D / 2 + 1 (a/r]) 1 



A(^) « 77T ^ T7 -— 7J Re[(2a/r ? ) 2 ^,(^)], (37) 



where we have defined the function 

OO r- 

g u {z) = T{u)T(D/2 - u) ^(<5i<5 2 ) n 2n 2 ^ D + Sj(n + \z- a 3 \/a) 2u ~ u . CW) 

n=l L J=l,2 



As it is seen from Eq. (|37p , for large separations of the plates the behavior of the interference 
term is qualitatively different for real and for imaginary values of the parameter v. For positive 
values of v , to leading order, we find 

For imaginary v, the asymptotic behavior at large distances is of the form 

A/ 2 ^ „ M2 D )|cos[2Hln(2aA ? ) + </>] 

> 2^TW(^p ' ^ 

with g v {z D ) = \g u (z D )\e ir f' . Hence, in this case the decay of the interference part is oscillatory. 
At a given spatial point, the dependence on the synchronous time coordinate has the form 
A(y? 2 ) ~ e -Dt/a cos [2\v\t/a + if}]. One gets similar oscillations for the single plate parts at large 
distances from the plate (see Ref. [9j). Note that, for a given a, the limit under consideration 
corresponds to the one for the late stages of the cosmological expansion. 



4 Vacuum expectation value of the energy- momentum tensor 



For the evaluation of the VEV of the energy-momentum tensor in the region between the plates, 
we use 



(0|T <fc |0) = lim did' k W(x,x') + [(£ - 1/4) g if Sjy - £V,V fc - CR ifc l (if 



(41) 



where Rik = Dgik/a 2 is the Ricci tensor for the dS spacetime. Taking advantage of the ex- 
pressions for the Wightman function and for the VEV of the field squared from the previous 



9 



sections, the renormalized VEVs for the diagonal components of the energy-momentum tensor 
can be expressed in the form (no summation over /) 



(Tl) = (T/) dS + {Tf) j + 



dyy 



l-D 



dx H(x, y) 



x [g(^x/rj, \z D - ajlx/^diy) + 2G lX 2 F v {y)} , 
where we have introduced the notations 
G (y) 



(42) 



[(yi/A)d 2 y -D(t + tD)yd y + D 2 Z 
+m 2 a 2 - y 2 + (1 - A^) x 2 ] F u (y), 



G D {y) 
Gi{y) 



C-- 4 )y 2 d 2 + 



y 2 — x 2 



e(2 -d) + 



D-l 



yd y -HD\F v {y) 



(43) 



G D (y) + 



D-l 



+ (l-4£)x 2 



and 



G D = 1, Gi = 4£-1, 1 = 0,1, 



F„(y), 1 = 1,..., D-l, 



D-l. 



(44) 



Note that, though not explicitly written, the functions Gi(y) with I = 0, . . . , D — 1 do depend on 
x as well. In Eq. (j4"2"j) . (Tj)dS is the corresponding renormalized VEV in dS spacetime without 
boundaries. For points away from the plates, renormalization is strictly necessary for this part 
only. Owing to the dS invariance of the Bunch-Davies vacuum, the part (Tj^ds is proportional 
to the metric tensor with a constant coefficient and has been well investigated in the literature 
[131 HI- F° r t ne P ar t induced by a single plate at z D = aj, one has [9] (no summation over I) 



A 



D 



D+l 



a 



dyy 



l-D 



o 



dx(x 2 -y 2 )( D -V' 2 



u 



2x\z D - aj \/ V 

x — Gi(y). 



Cj(x/rj) 

The last term on the rhs of Eq. (|42p is induced by the presence of the second plate. 
For the non-zero off-diagonal component, we have 



(45) 



®5 



sen z 



D 



A D 



2a D ~ 



roo 

T / dyy 1 - D G D(y) 
Jo 



where the part corresponding to the geometry of a single plate is given by 



x / dxxH(x,y) Cj {xlr})e 2x \ zD - a ^ -e- 2x \ zD - a ^/ Cj {x/ri) 



(46) 



sgn(z 



D 



dyy l - D G QD (y) 



x / dxx(x 2 -y 2 )^' 2 



e -2x\z u - aj \/r) 
Cj(x/f]) 



In these formulae we have defined the function 

G 0D {y) = m-l)yd y + Ai]F v {y). 



(47) 



(48) 



The off-diagonal component (|46p corresponds to the energy flux along the direction perpendicular 
to the plates. This type of the energy flux also appears in the geometry of a cosmic string on 
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backgrounds of Friedmann-Robertson- Walker and dS spacetimes [TS]. Depending on the values 
of the coefficients in the boundary conditions and of the field mass this flux can be positive or 
negative. As an additional check of the expressions for the energy-momentum tensor, it can be 
seen that the ones for the single plate contribution and for the second plate induced part fulfill 
the trace relation 

(T/> = [d(£ - ^)V / V Z + m 2 ] <^>. (49) 

The boundary induced part in the VEV of the energy-momentum tensor is traceless for a con- 
formally coupled massless scalar. The trace anomaly is contained in the boundary-free part 
only. 

For a conformally coupled massless scalar field (£, = m = 0) the single plate part in the 
VEV of the energy-momentum tensor vanishes and one finds 

= <^- (4j )X°W2+l) 

f°° x D 
x / dx — — — -. (50) 

As in the case of the field squared, the boundary induced part in this formula could have 
been obtained from the corresponding result for the Casimir effect in Minkowski spacetime, 
by using the fact that the two problems are conformally related. The electromagnetic field 
in D = 3 is conformally invariant and the Casimir problem with two perfectly conducting 
parallel plates is reduced to the corresponding problem with two scalar modes with Dirichlet 
and Neumann BCs. In this case, the single plate parts vanish and for the interference part 
we have A(7] fc ) = — (7r 2 /720)(aa/?7)~ 4 diag(l, 1, 1, — 3). In the case D > 3, the electromagnetic 
field is not conformally invariant and we expect that the corresponding VEV will depend on the 
distance from the plates. However, this case requires further consideration. 

The VEVs for the components of the energy-momentum tensor can be written in the more 
symmetric form 

(T l k ) = (T}) dS +J2(T l k ) j + ^(Tl), (51) 

3=1,2 

where for the interference terms we have (no summation over I) 

A(T/> = J dyy l - D J dxH(x,y) 

x [h(x/ V , z D )G l (y) + 2G lX 2 F u {y)] , (52) 

A pOO /"CO 

A(T D ) = ^t/ o dyy l - D J dxxH(x,y) 

e -2x\z D - aj \/r, 

x sgn(z - aj ) G 0D {y). (53) 

i=i,2 cjWn) 

Note that when the coefficients in the BCs are the same, (3\ = 02, the energy flux vanishes at the 
point z D = (a\ + 02) /2. Of course, this is a direct consequence of the symmetry of the problem. 

In the special cases of Dirichlet and Neumann BCs, expressions similar to (|33p can be ob- 
tained for the interference terms in the VEVs of the energy-momentum tensor. Here we give the 
corresponding formulas for the interference parts in the normal stress and in the off-diagonal 
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component: 



A ( T d) = (27r) D/ 2+ i E I dy^y[{D-l)f D/2 {z) + f D/2 _ 1 {z)]F v {y) 



<5( J ) D 
fD/2-i( z ) - — 2^ fD/2-i{z + %y\z ~ aj\/v) 
i=i,2 



^} , (54) 

£* J z=2yna/r) 



xz fD/2(z)\z=2y(an+\z D ~ aj \)/T,GoD(y), (55) 



with J=D,N and being the function f^,(z) defined after the formula (|33p . As for the case of the 
field squared, the coordinate dependent parts have opposite signs for Dirichlet and Neumann 
boundary conditions, respectively. 

We now wish to examine the behavior of the VEV of the energy-momentum tensor in the 
asymptotic regimes of small and of large separations between the plates. At small separation, 
assuming that a/i] <C 1, we introduce in (|52p . (|53|) new integration variables u = ax/rj and 
v = ay/rj. The arguments of the functions Gi{y) and F v (y) are large and we can use the 
corresponding asymptotic formulae for the modified Bessel functions. In particular, one has 
Fu{y) ~ y D ~ l and 

G (y) « [(1 - 40 x 2 - y 2 ] y D -\ G D (y) « D& - 0y D ~\ 

Gi(y) « G (y) + Dy *~f y D -\l = l,...,D-l. (56) 

For the further evaluation of the integrals, we use Eq. (|35p with v = ±1/2. As a result, to 
leading order, one has (no summation over I) 

a{t!) k r^m^+cj 



T(D/2)(a/r]) D + 1 J c 1 (x)c 2 {x)e 2ax - 1 ' 

D „ 2(4n)- D / 2 D(Z - gp) v sgn^-q^x^^e-^l^-^l 
/ ~ T(D/2)a(a/r,)0 ^J X Cl (x)c 2 (x)e^ - 1 c,(x) ' 



where 



Bl = "2 (? - Cd) , B D = 0, = 4£ - 1, C D = 1, (58) 

with I = 0, 1, . . . , D — 1. Note that A(Tq ) ) ~ (a/rj)A(T}) and that, for a conformally coupled 
field, the leading terms in the diagonal components are homogeneous. For small separation, 
to leading order, the energy density is equal to the stresses along the directions parallel to the 
plates. The vacuum stress normal to the plates, A(T^), is positive for Dirichlet and Neumann 
BCs, and is negative for Dirichlet BC on one plate and Neumann on the other. 

Now let us discuss the asymptotics at large distances between the plates, a/r] 3> 1, when the 
curvature effects are essential. The corresponding asymptotic formulae for the the interference 
parts in the VEV of the energy-momentum tensor are found in a way similar to that already 
described for the VEV of the field squared in Sect. OH For positive values of the parameter v, 
the leading terms have the form (no summation over I) 

A<r/> . f«*- D - x 9^ D ) 



47r D/2+i {2a/r]) D - 2v ' 
47r D/2+i (2a/r]) D - 21 ' - J ' 



A<*iP> - aJ^J:±> > (59) 
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where the function g u {z) is defined by Eq. (|3"5|) and 

$Hz) = [(4$-l)(D-2v)+4$]r(D/2-i, + l) 

n=l ,7=1,2 v 1 5 u ' 

The coefficients /; for separate components of the energy-momentum tensor are defined as 

/, = (2v/D)f = -2v [£ + (£- 1/4) (£> - 2i/)] , (61) 

with I = 1,...,D. Note that for minimally and conformally coupled massless fields /o = 
fl = and thus the leading terms vanish. As we see, in the limit under consideration the 
vacuum stresses are isotropic and |A(Tq )| <C |A(T/)|. The corresponding equation of state is 
of barotropic type: A(T 1 1 ) w • • • w A(Tg) « (2v/D)A(T$). A similar relat ion holds for the 
single plate parts at large distances [9]. As a consequence, the equation of state parameter for 
the total boundary-induced part is equal to —2u/D, and it is negative. 
For imaginary v the leading terms at large distances are given by 

2^/^+1 (2a/r/) 

A(T ° D) W 2,Ci^Ji COS [ 2 '^ ln ^ 2Q /^ + ^ (62) 

In these formulae // = |/z|e 1 ^' and 5^(2"°) = |<?^(2' D )|e* l ^° . In this case, the damping of 
the interference parts as functions of the proper distance is oscillatory. In terms of the syn- 
chronous time coordinate, at a given spatial point we have damping oscillations in accordance 
with A(T/) ~ e Dt / a cos(2\v\t/a + if>i). From (]6ip it follows that the oscillations in the energy 
density and in the vacuum stresses are shifted in phase by ir/2. 

In Fig. [U for the case of a D = 3 conformally coupled scalar field with Dirichlet BCs on both 
plates, we have plotted the energy flux as a function of z D /r), for given values of a/rj = 1 and 
ma = 1/4, 1 (left plot), and as a function of ma for given values of ajr\ = 1 and z D /rj = 0.3 
(right plot). For Neumann BCs the flux has the opposite sign. 



5 Casimir force 



In this section we consider the vacuum forces acting on the plates. The vacuum force acting per 
unit surface of the plate at z D = a,j is determined by the ^-component of the vacuum energy- 
momentum tensor evaluated at this point. For the region between the plates, the corresponding 
effective pressures can be written as a sum of two terms, namely 

pV) =p? ) +pj£ t) , J = 1,2. (63) 

The first term on the rhs is the pressure for a single plate at when the second plate 

is absent. This term is divergent due to the surface divergences in the subtracted vacuum 
expectation values and needs additional renormalization. The second term on the rhs of Eq. (|63p 
is the pressure induced by the presence of the second plate, and can be termed as an interaction 
force. This contribution is finite for all nonzero distances between the plates. In the regions 
z D < a\ and z D > a,2 we have 

p (j) = p (j\ 

As a result, the contributions to the vacuum force 
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z M ma 



Figure 1: Energy flux in the region between the plates as a function of z D /rj for a/rj = 1, 
ma = 1/4,1 (left plot), and as a function of ma for a/i] = 1, z /rj = 0.3 (right plot), for a 
D = 3 conformally coupled scalar field with Dirichlet boundary condition. On the left plot the 
figures near the curves are the values of the parameter ma, and the scaling factor n is equal to 
2 and 3, for ma = 1 and ma = 1/4, respectively. 

coming from the term ' are the same from the left and from the right sides of the plate, so 
that there is no net contribution to the effective force. 

The interaction force on the plate at zP = aj is obtained from the last term on the rhs of 
Eq. (|4*2j) for (Tp) (with minus sign) substituting z D = a,j, that is 

= f°° dyy 1 ' f°° dxx 2 H{x,y) 

a Jo Jy 

(64) 



(Pjx/vY - 1 

Depending on the values of the coefficients in the boundary conditions, the effective pressures 
(|64"|) can be either positive or negative, leading to repulsive or to attractive forces, respectively. 
For j3\ /?2 the Casimir forces acting on the left and on the right plates are different. In the 
Appendix we show that, as it must be, in the limit a — > oo the corresponding result for the 
geometry of two parallel plates in Minkowski spacetime is obtained. 

The general formula is further simplified for the special cases of Dirichlet and of Neumann 
boundary conditions. For Dirichlet BCs on both plates, we find 

A —D—l 00 roo 

= - (27r )D/2+i E/ d yy p M i( D -l)fD/2&nay/n) + f D /2-i(2nay/ V )] . (65) 

For ^ v < 1 the integrand in this formula is positive and is negative, yielding an attractive 
force for all separations. In the case of Neumann BCs the corresponding expression is 

{N ) (D) 8a- - 1 ^ (°° G D (y) 
P(int) = P(int) - ( 27r )D/2 + i 2-, J Q dy f D /2 ^(2nay / V ) , (66) 

where the function Go(y) is defined in Eq. (j4"3~|) . 
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Let us now investigate the asymptotic behavior of the vacuum forces. This can be done in the 
way similar to that already described in the case of the VEV for the energy-momentum tensor. 
In the limit of small proper distances between the plates, a/ij <C 1, the main contribution to the 
integral over y comes from large values of y, y ~ Tj/a. By using the corresponding asymptotics 
for the modified Bessel functions, to leading order we find 

(i) 2(r ? /a) g+1 f°° x D 

P(int)~ (47r) D / 2 r(L»/2) J Cl (x)c 2 (x)e 2ax -I [ ' 

If, in addition, \(3j\/a 3> 1, one has 

W ^ OT((g + l)/2)C R (g + l) , . 

^M)~ ( 4vr )(D+i)/2 (aa/^D+i ' 

and the corresponding force is attractive. In (|68p . Cr( x ) is the Riemann zeta function. The same 
result is obtained for Dirichlet BCs on both plates. In the case of Dirichlet BC on one plate and 
non-Dirichlet one on the other, the leading term is obtained from (j68H with an additional factor 
(2~ D — 1). In this case the vacuum force at small distance is repulsive. 

Now we consider the large distance asymptotics, a/rj 3> 1. The cases of real and imaginary 
v must be studied separately. For positive values of v, one has 

(i) 2a~ D ~^ 

P(int) ~ 7r D/2+l( 2o /^)0-2i/+2' 1 < 00 ' 

(i) or"" 1 ^ 

~ " 7 r^/ 2 + 1 (2a/7 ? )^ 2 -' ^ = °° ' (69) 

where we have introduced the notations 

g® = [(D + l)/2-u}T(D/2-u + l)T(u) 

x[i-2w j/V ) 2 f D }f:^%, ( 7 °) 

n=l 

^Ci) = r(D/2-,)rH/ D f; (5l52)n 



n=l 



n D-2v 



and /d is defined by Eq. (j6T|) . Note that (3j = oo corresponds to Neumann BC. In the case of 
non-Neumann BCs we have assumed that \f3j\/a <C 1. As it is seen from (|69h . when fjj ^ 0, 
at large distances the ratio of the Casimir forces acting on the plate with Neumann and non- 
Neumann BCs is of the order (a/ry) 2 . Note that in neither of these cases does the force depend 
on the specific value of Robin coefficient in the BC on the second plate. For Dirichlet BC on 
the plate at z D = aj (f3j = 0), at large separations the Casimir force acting on that plate 
is repulsive (attractive) for Neumann (non-Neumann) BCs on the other plate. The nature of 
the force acting on the plate with Neumann BC depends on the sign of fo and can be either 
repulsive or attractive, in function of the curvature coupling parameter and of the field mass. 
For minimally and conformally coupled massive scalar fields one has fjj = v(D/2 — v) and 
fo = v(l/2 — v)/D, respectively, and this parameter is positive. The corresponding force is 
attractive (repulsive) for Neumann (non-Neumann) BC on the second plate. Note that for the 
geometry of parallel plates in the Minkowski bulk the Casimir forces at large distances are 
repulsive for Neumann BC on one plate and for non-Neumann BC on the other plate. For all 
other cases of BCs the forces are attractive. 
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At large separations between the plates and for imaginary v, the leading order terms in the 
corresponding asymptotic expansions take the form 



0) 

P(int) 



JJ) 
P(int) 



4a- D " 1 |^ ) | 



• cos[2|^| ln(2a/??) + 4>{j)], < °°, 



TT D / 2 + 1 (2a/r ] ) D + 2 
2aT D - 1 \a m \ 



oo, 



(71) 



where we have defined gl~ 



'W) and 



In such case the decay of the 



vacuum forces is oscillatory. In terms of the synchronous time coordinate one gets the behavior 
p^l t) ~ exp[(D + 2 - 25 0>1/ p.)t/a] cos[2\u\t/a + ^J. 

In Figs. [2] and [3] we have plotted the Casimir force for a D = 3 scalar field, conformally 
and minimally coupled to gravity, respectively, as a function of the proper distance between the 
plates, measured in units of the dS curvature scale, a. The left (resp. right) plots are for Dirichlet 
(resp. Neumann) BCs on both plates. The figures near the curves correspond to the values of the 
parameter ma. Values are taken in a way so to have both possibilities, with positive and purely 
imaginary values of the parameter v, with corresponding monotonic and oscillatory behavior of 
the forces at large distances, respectively. Note also, in particular, the plots corresponding to a 
massless field in the two cases, corresponding to a photon-like contribution. 
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Figure 2: Interaction forces between the plates for a D = 3 conformally coupled scalar field with 
Dirichlet (left plot) and Neumann (right plot) BCs. The figures near the curves are the values 
of the parameter am. Note in particular the plots corresponding to a massless field in the two 
cases. 

In Fig. H] the dependence of the Casimir force on the parameter ma for a given separation 
corresponding to a/r] = 3 is depicted. Conformally coupled scalar fields with Dirichlet and 
Neumann BCs, respectively, are considered. 

From the discussion given above it follows that for the proper distances between the plates 
larger than the curvature radius of the dS spacetime, aa/r] > a, the gravitational field essentially 
changes the behavior of the Casimir forces compared with the case of the plates in Minkowski 
spacetime. In particular, the forces may become repulsive at large separations between the 
plates. In particular, for real values v and for Neumann BC on both plates, Casimir forces are 
repulsive at large separations, in the range of parameters for which fo < [see Eqs. ([69]) . ([70]) ]. 
Recall that, for the geometry of parallel plates on the background of Minkowski spacetime, 
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Figure 3: Same as in Fig. [2] for a minimally coupled scalar field. 




Figure 4: Interaction force between the plates for a/r] = 3 as a function of the field mass, for a 
D = 3 conformally coupled scalar field with Dirichlet and Neumann BCs. 
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the only case with repulsive Casimir forces at large distances corresponds to Neumann BC on 
one plate and non-Neumann BC on the other. A remarkable feature of the influence of the 
gravitational field is the oscillatory behavior of the Casimir forces at large distances, which 
appears in the case of imaginary v. In this case, the values of the plate distance yielding zero 
Casimir force correspond to equilibrium positions. Among them, the positions with negative 
derivative of the force with respect to the distance are locally stable. 

In the discussion above we have considered the expectation values assuming that the field is 
prepared in the Bunch-Davies vacuum state. This corresponds to the effect of vacuum polar- 
ization by boundary conditions. In states containing particles the expectation values of physical 
observables will receive additional contributions. For example, the expectation value of the field 
squared has the form (</? 2 ) = (</J 2 )bd + 2 J2i 1 " n o-(Pa{x) l Pa(x), where l n a is the number of particles 
with the set of quantum numbers a and (</? 2 )bd is the expectation value in the Bunch-Davies 
vacuum state. On the base of this relation, we can consider the effects of the boundaries at 
finite temperature assuming that the field is in thermodynamical equilibrium. In this case, the 
VEVs are changed by the thermodynamical expectation values. However, it should be noted 
that, owing to the time dependence of the background spacetime, we can talk about thermo- 
dynamical equilibrium in the adiabatic limit only. In dS spacetime this corresponds to the 
conditions r/K 3> ma and r/K 3> 1, for the modes with wave number K (see also the discussion 
in |14|). By taking into account that, at temperature T, the main contribution to the thermo- 
dynamic expectation values comes from the region K < T, we obtain the conditions T c ^$> 1/a 
and T c m, where T c = r/T/a is the comoving temperature. The dominant contribution to 
the boundary induced expectation values comes from the fluctuations with K < 1/a. Combin- 
ing this with the estimates given above, we conclude that the adiabatic approximation for the 
boundary induced expectation values corresponds to the limit a/rj <C 1. As has been shown 
above, in this limit the leading terms in the VEVs induced by the plates coincide with the 
corresponding quantities for the geometry of parallel plates in Minkowski spacetime. The same 
is true for the thermal corrections. 

6 Conclusion 

Amongst the most interesting topics in the investigation of the Casimir effect is the dependence 
of the characteristics of the vacuum fluctuations on the background geometry. In the present 
paper we have considered the classical geometry of two parallel plates on the background of 
dS spacetime for a scalar field with Robin boundary conditions on the plates. The general 
case has been investigated when the constants in the Robin boundary conditions are different 
for the two separate plates. In the region between the plates, the Wightman function has 
been obtained and displayed under the form of a mode sum involving series over zeros of the 
function defined by Eq. ([8]). For the summation of this series we have made use of expression 
()12p . This has allowed us to extract, from the Wightman function, the part coming from a 
single plate, and to express the additional part in terms of integrals, which are exponentially 
convergent in the coincidence limit. The single plate contribution was investigated previously, in 
Ref. [9]. The contribution induced by the second boundary has been presented in two alternative 
forms, as given by Eqs. (]14p and (I22D . By using the expression of the Wightman function, we 
have evaluated the VEVs of the field squared and of the energy-momentum tensor, in the region 
between the plates. These VEVs are decomposed into a boundary-free dS, a single plate-induced 
and an interference contribution, respectively. The last one, for the cases of the field squared 
and energy-momentum tensor, is given by Eqs. (i30j) and (i52j) . (i53j) . respectively. The vacuum 
energy-momentum tensor is non-diagonal, with the off-diagonal component corresponding to 
the energy flux along the direction normal to the plates. In the case of a conformally coupled 
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massless field, the total single plate contribution to the VEV of the energy-momentum tensor 
vanishes and the interference part is obtained from the corresponding result for the Minkowski 
bulk, by standard conformal transformation. 

Various limiting cases have been studied. In the limit of small distances between the plates 
the interference part in the VEV of the field squared coincides, to leading order, with the corre- 
sponding quantity for a conformally coupled massless field, and is given by Eq. (|32p . The leading 
terms of the interference parts of the VEV for the energy momentum tensor are given by expres- 
sions (|57p . For a conformally coupled scalar field, the leading term of the off-diagonal component 
vanishes, and the leading terms of the diagonal components are homogeneous. In the opposite 
asymptotic limit of large separations between the plates, the behaviors of the interference parts 
crucially depend on the value of the parameter u, defined by Eq. ((7|). For positive values of this 
parameter, the leading terms of the corresponding asymptotic expansions are given by Eqs. (J39]) 
and (|59p . for the field squared and the energy-momentum tensor, respectively. The interference 
contributions for the field squared and the diagonal components of the energy-momentum ten- 
sor decay as l/(a/r]) D ~ 2u , whereas the off-diagonal component decays like l/(a/r)) D ~ 2u+1 . To 
leading order, the vacuum stresses are isotropic and the boundary induced VEV of the energy- 
momentum tensor corresponds to a gravitational source of barotropic type, with equation of state 
parameter equal to — 2v/D. At large separations between the plates and for imaginary values of 
the parameter v, the asymptotic behavior of the interference parts for the field squared and for 
the energy- momentum tensor are given by Eqs. (140p and ([62]), respectively. The corresponding 
behavior is damping oscillatory and the VEVs decay as (a/r])~ D cos[2|z^| In (2a/rj) + ift], for the 
field squared and the diagonal components of the energy-momentum tensor. For the off-diagonal 
component the amplitude decays as (a/t])^ D ^ 1 . 

The vacuum forces acting on the plates are determined by the ^-component of the stress. 
They have been studied in Sect. EJ The normal stresses on the plates are presented as sums 
of single plate and interaction contributions. The contributions to the vacuum force coming 
from the single plate terms are the same from the left and from the right sides of the plate 
and thus give no contribution to the effective force. The interaction forces per unit surface are 
determined by formula (|64p . This expression is further simplified in the special cases of Dirichlet 
and Neumann BCs, yielding Eqs. (j65|) and ([66]) . respectively. For small distances between the 
plates, to leading order the vacuum forces are given by Eq. (j67l) . If, in addition, \(3j\/a S> 1, 
the vacuum forces are attractive at small distances, except for the case of Dirichlet BC on one 
plate and non-Dirichlet on the other, in which case the force turns out to be repulsive. At 
large distances between the plates and for positive values of v, the force acting on the plate 
decays monotonically as l/(2a/r]) D ~ 2u+2 , for non-Neumann BCs, and as l/(2a/rj) D ~ 2u , in the 
case of Neumann BCs [see Eqs. (]69p ]. For imaginary values of v the behavior of the vacuum 
forces is damping oscillatory, in the leading order described by Eqs. (|7ip . Having in mind that 
spectral properties of spin 2, 1, operators for dS spacetime are known, the current study can 
be now extended to the calculation of the Casimir force due to quantum gravity (for the one- 
loop effective action of arbitrary quantum gravity in dS spacetime see Ref. |21|). This will be 
considered elsewhere. 

From the analysis carried out above, it follows that the curvature of the background space- 
time decisively influences the behavior of boundary induced VEVs at distances larger than the 
curvature scale. As we have seen, when the background is dS spacetime the decay of the VEVs at 
large separations between the plates is power-law (monotonical or oscillating), independently of 
the field mass. This is quite remarkable and clearly in contrast with the corresponding features 
of the same problem in a Minkowski bulk. To wit, the interaction forces between two parallel 
plates in Minkowski spacetime at large distances decay as l/a D+1 for massless fields and these 
forces are exponentially suppressed for massive fields by a factor of e~ 2ma . For the geometry of 
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two parallel plates, in AdS spacetime the decay of the vacuum forces at large separations is also 
exponential (see Ref. [3]), with the suppression factor being determined by the AdS curvature 
scale. In a way very much similar to the procedure described in Ref. [22J (see also Ref. |23] 
for finite temperature effects), we are able to treat here the more general case of dS spacetime 
with compact internal subspaces and piston-like geometries. Note that this calculation can be 
extended to a self-interacting scalar field theory too, in which case mass becomes an effective 
mass, proportional to the background scalar. In this way, our results and procedures here can 
be used to the study of curvature-induced phase transitions of the in-in effective potential in the 
same way as it was proposed for the out- in effective potential in Ref. [24J. 
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A Minkowski spacetime limit 

In this appendix we will explicitly show the limiting transition of the situations considered above 
to the geometry of two parallel Robin plates in Minkowski spacetime, for the vacuum interaction 
forces. In this limit a — > oo and the modulus of the order of the modified Bessel functions is 
large, v ~ ima. In addition, we have r\ — > a. Changing the integration variables to x = ur), 
y = vrj, we see that the arguments of the modified Bessel functions are large, too. Hence, 
we make use of the uniform asymptotic expansions for these functions for imaginary values of 
the order with large modulus. The leading terms in these expansions have the form (see, for 
example, [20J) 

K^nz) ~ TV^T^/ 2 cos[/i/(z///) - vr/4], 

2g/i7r/2 

sin[/i/(z/^)-7r/4], (72) 



for z < 1 and 



for z > 1. The functions in these formulas are defined as 



1 + VT 



2" 



/(*) = M ; 1-vT 

g(z) = —arcsec z + \J z 1 — 1. (74) 
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From Eqs. ((72} and {7J} it follows that 



I v {y)K u (y) ~ cos [2maf(y/ma) ], y < ma, 

ma 

1 y D 

I u {y)K u {y) ~ =, y > ma. (75) 

ma v(y/ ma ) ~ 1 

The main contribution to the force (|64p comes from the region (ma, oo) of the integration over 
y; to leading order, we find 

W „ 2(4vr)-^ r« M 2 (M 2_ m 2 )g /2-l 

p M) ~ rp/2) i m (giu-i)(fau-i) e2att _ 1 • W 

(/3i«+l)(/3 2 ii+l) 

This result coincides with the corresponding formula for parallel plates in the Minkowski bulk. 
In a similar way, the limiting transition for the other quantities can also be checked. 
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